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Abstract 

X 

C^ I The electroweak process pp — > £^1/77 is calculated at tree level, including 

finite W width effects. In order to obtain a gauge invariant amplitude, the 
imaginary parts of WW'j triangle graphs and WW'yy box diagrams have to 
be included, in addition to resumming the imaginary contributions to the 
W vacuum polarization. We demonstrate the existence of a radiation ampli- 
tude zero in pp -^ W^^j -^ £^1^77, and discuss how it may be observed in 
correlations of the 77 and lepton rapidities at the Fermilab Tevatron. 



I. INTRODUCTION 

Radiative W production and decay at hadron colliders is an important testing ground 
for the Standard Model (SM). The simplest process, qq' -^ i'^uj, allows the measurement 
of the WW'y three gauge boson coupling at large photon transverse momenta HTj-Q. In 
addition, this process is of special interest due to the presence of a zero in the amplitude of 
the parton level process qq' — > W'~f ||T|J§,§]- At small transverse momenta of the photon or 
when the photon is emitted collinearly to the final state charged lepton, this process needs 
to be fully understood when trying to extract a precise value of the W boson mass from 
Tevatron data. Approximately 24% (13%) of all W ^ eu {W -^ fiu) events contain a photon 
with a transverse momentum (pj,) larger than 100 MeV 0J^, the approximate threshold of 
the electromagnetic calorimeter of the CDF and D0 detectors. Radiative W decay events 
shift the W mass by about 65 MeV in the electron, and by approximately 170 MeV in the 
muon channel P,p!0|. 

For similar reasons, the process qq' — > i^wy'y is interesting. At large photon transverse 



momenta, W'j'j production is sensitive to the structure of the WW'j'j quartic coupling [jTT 
Furthermore, as a consequence of a general theorem [^ one expects a radiation zero in the 
SM qq' — > W'-ff amplitude. Two photon radiation is expected to have a non-negligible 
effect on the W mass extracted from future high precision Tevatron data because approxi- 
mately 0.8% of all W -^ fiu events are expected to contain two well separated photons with 



Pj: > 100 MeV [T^. Finally, W'j'j production is an irreducible background to associated 
production oi a W and a Higgs boson in hadronic collisions, if the Higgs boson decays into 
two photons [|l^. 

In pp collisions at a center of mass energy of 2 TeV, the total cross section for W^jj 
production is approximately 4.6 fb, when only considering leptonic decays, W ^ iu {i = 
e, /i), and pj^ > 10 GeV and |?7^| < 2.5 {rj being the pseudorapidity) |1^. Upgrades of 
the Tevatron accelerator complex (TeV33), beyond the Main Injector project, could yield 
an overall integrated luminosity of (9(30 fb"^) |[T^, making a study of PF77 production a 



realistic goal in the TeV33 era. The hadronic decay modes of the W will be difficult to 
observe due to the QCD JJ77 background |T^. We therefore concentrate on the leptonic 
decays of the W boson, and calculate the helicity amplitudes for the complete process 

qq' -^ i^u-fy, (1) 

including Feynman diagrams where one or both photons are emitted from the final state 
charged lepton line. In a realistic simulation, these diagrams, together with finite W width 
effects, need to be taken into account. 

When including finite W width effects, some care is needed to preserve gauge invariance. 
Replacing the W propagator, l/(g^ — m^), by a Breit-Wigner form, 1/ {q"^ — m'^ + im]yT]y) , 
will disturb the gauge cancellations between the individual Feynman graphs and thus lead 
to an amplitude which is not electromagnetically gauge invariant. In addition, a constant 
imaginary part in the inverse propagator is ad hoc: it results from fermion loop contribu- 
tions to the W vacuum polarization and the imaginary part should vanish for space-like 
momentum transfers. In Ref. [^ it was demonstrated how this problem is solved for W'y 
production by including the imaginary part of WW'y vertex corrections in addition to the 
resummation of the W vacuum polarization contributions. Here, we generalize the result 
of Ref. P] to W'jj production, and show that a gauge invariant amplitude for the process 
qq' — > i'^u'-f'-f is obtained by also including the imaginary part of the WW'-ff box correc- 
tions. Extending the argument, a gauge invariant amplitude for qq' ^ W + w-f, n > 2, 
can be obtained by implementing the corrected WW'j and WW'j'y vertices together with 
the resummed W vacuum polarization contributions. No higher WWri'y vertex functions 
need to be considered. Our analysis of gauge invariance for W'j'j production is described in 
Sec. II. 

The existence of a radiation zero in the process qq' — > W'y has been well known for more 
than fifteen years |^,p . All SM helicity amplitudes for the process qq' — > W'y vanish at 

cos 9^ = cos 9*Qw = " _ "' , (2) 



where 6^ is the angle between the W and the incoming quark q, in the parton center of 
mass frame. A theorem then predicts that the process qq' ^ W + wf, n > 1, exhibits a 
radiation zero for the same scattering angle cos 6'qvi/, if the n photons are coUinear. In Sec. Ill, 
we numerically demonstrate the existence of this radiation zero in Wj'j production. 

In practice radiation zeros in hadronic collisions are difficult to observe. In the W'y 
case, the ambiguity in reconstructing the parton center of mass frame and in identifying the 
quark momentum direction represents a major complication in the extraction of the cos^^ 
distribution f^. Higher order QCD corrections [IT^jH] and finite W width effects, together 
with photon radiation from the final state lepton line, transform the zero to a dip [jl9 



Finite detector resolution effects further dilute the radiation zero. The twofold ambiguity 
in reconstructing the parton center of mass frame originates from the nonobservation of the 
neutrino arising from W decay. Identifying the missing transverse energy with the transverse 
momentum of the neutrino, the unobservable longitudinal neutrino momentum, Pl(i^), and 
thus the parton center of mass frame, can be reconstructed by imposing the constraint 
that the neutrino and charged lepton four momenta combine to form the W rest mass [ED 



The resulting quadratic equation, in general, has two solutions. Finally, determining the 
cos^^ distribution requires measurement of the missing transverse energy in the event. In 
future Tevatron runs, one expects up to ten interactions per bunch crossing [^. Multiple 



interactions per crossing significantly worsen the missing transverse energy resolution, and 
thus tend to wash out the dip caused by the radiation zero. 

For W'j'y production, the same problems arise. In addition, it is very difficult to experi- 
mentally separate two coUinear photons, and, thus, to distinguish the W'yy signal from W'y 
events and from the W+ jets background, where one of the jets fiuctuates into a 7r° which 
decays into two almost coUinear photons. One therefore has to search for a signal of the 
radiation zero which survives an explicit photon-photon separation requirement. 



In Ref. [^] it was found that lepton-photon rapidity correlations offer the best chance 
to observe the radiation zero in W'y production. The distribution of the rapidity difference 
A|/(7,£) = iiy — iii clearly displays the SM radiation zero. It does not require knowledge 



of the longitudinal momentum of the neutrino, and so automatically avoids the problems 
described above. In Sec. Ill we show that the concept of rapidity correlations as a tool to 
search for radiation zeros can be generalized to the W'j'j case. The A?/(77,£) distribution 
with cos6'-y^ > 0, where 6'^^ is the opening angle between the two photons in the laboratory 
system, clearly displays the SM radiation zero even when one requires two well separated 
photons, provided that cuts are imposed which reduce the background from radiative W 
decays. 

It is sometimes useful to compare distributions for W'j'y and ^'77 production. Simulta- 
neously with the calculation of the process qq' —>■ i^wyj, we therefore also present results 
for qq -^ i'^i~'y'-f production in Sec. III. Section IV contains some concluding remarks. 

II. FINITE WIDTH EFFECTS AND GAUGE INVARIANCE IN ^^77 

PRODUCTION 

At the parton level, the reaction pp — > i^w-f-f proceeds via the Feynman diagrams shown 
in Fig. |1|. Besides the diagrams for W'j'j production, graphs describing W'j production 
followed hj W ^ iwj contribute, and also the radiative decay W — ^ iwj'j. When finite W 
width effects are included, the three reactions can no longer be distinguished, and the full set 
of Feynman diagrams must be taken into account. To calculate the helicity amplitudes for 
the process qq' — > i^wj'j we have used the framework of Refs. p2| and pSf. The result was 



then compared numerically with the amplitudes obtained using the MADGRAPH/HELAS 
program ||2^ , |25[| which generates helicity amplitudes automatically. All quarks and leptons 
were assumed to be massless in our numerical calculations. 

A naive implementation of finite W width effects, by replacing the W propagator by a 
Breit-Wigner form with momentum dependent width, leads to a violation of electromagnetic 



gauge invariance [^ and the resulting cross sections cannot be trusted. One encounters the 



same problem when gauge invariance is restored in an ad hoc manner [HlpTl- Finite width 



effects are included in a tree level calculation by resumming the imaginary part of the W 
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FIG. 1. Feynman diagrams for the process qq' -^ ii/j'j. Permutations of the final state 
photons and the W boson are not shown explicitly. 



vacuum polarization. Gauge boson loops (W^7 and WZ) only contribute above the W- 
mass pole and are suppressed by threshold factors ||2^. They can safely be neglected at 
the desired level of accuracy and only fermion loops need to be considered. Neglecting the 
fermion masses in the loops, the transverse part of the W vacuum polarization receives an 
imaginary contribution 

i,nK(r)-Ei;,'--,'^_-,h., (3) 

while the imaginary part of the longitudinal piece vanishes. In the unitary gauge the W 
propagator is thus given by 

^ q'^ - rri^r + ilmll'^{q'^) \ (f J ni^ - ilmll^r{q^) (f 



^m-_:l^(i + ,^^) . (4) 



q^ - m^ + iq'^'jw V m^ j 

A gauge invariant expression for the amplitude of the process gg' — > £^z/77 is obtained 
by attaching the final state photons in all possible ways to all charged particle propagators 
in the Feynman graphs. To be specific, we shall concentrate on the £"1/77 final state in the 
following. In addition to radiation off the external fermion lines and radiation off the W 
propagators, the photons must be attached to the charged fermions inside the W vacuum 
polarization loops, leading to the fermion triangle and box graphs of Figs. |^ and ^ Since we 
are only keeping the imaginary part of 11^ (g^), consistency requires including the imaginary 
parts of the triangle and box graphs only. These imaginary parts are obtained by cutting 
the triangle and box graphs into on-shell intermediate states in all possible ways, as shown 
in the figures. 

For the momentum flow displayed in Figs. ^ and ^, the tree level WW'-f and WW'j'y 
vertices are given by the familiar expressions 



- zeTf^ = -ze {{q, + q^fg'^^ - {q, + kfg^^ + {k - q^Tg^^) , (5) 

-^eTf^'' = -^^ (2g-Pg^^ - g^^g^^ - g^^ g^^) . (6) 
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FIG. 2. The effective WW'j vertex as needed in the tree level calculation of ^1/77 production. 



For the triangle graphs of Fig. ^the niomentum configuration, ql > q^-, k'^ = 0, is the same 
as the one encountered in the qq' — » £^i"y case. Neglecting fermion masses, the full WW'j 
vertex is then given by ^ 



pa/3^ = r;^^^ 1 + E 



/ 



487r 



Tf^{l + t- 



mw 



Tf^ (1 + Z7H/) 



(7) 



Non-zero fermion masses, mf > 0, introduce corrections to Eq. (0) and generate axial 
vector contributions to the WW'-f vertex which are proportional to rn'j/ql and m?Jq2 p9| . 
They can be neglected at the desired level of accuracy for the lepton and the light quark 
doublets. Top-bottom loops do not contribute to the imaginary part of the WW'j vertex 
below threshold, i.e. for qf < {nit + nihY . In the imaginary part of the WW''yy vertex they 
are either absent for qf < {rrit + rn^Y or are suppressed by powers of the top quark mass. 
These massive loops are not needed for the restoration of electromagnetic gauge invariance 
and can be neglected close to the W pole. 

An evaluation of the 24 cut box diagrams of Fig. ^ yields a result similar to that of 
Eq- (0) 10 • For vanishing fermion masses, each fermion doublet /, irrespective of its 
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FIG. 3. The Feynman graphs contributing to the effective WWjj vertex of Eq. (^. Only 
four of the 24 cut box diagrams which contribute to the imaginary part of the one-loop WW'j'j 
vertex are shown. 



hypercharge, adds z(f7^/487r)ro '^ to the tree level WW'yy vertex Fg ''. In the phase space 
region qf > ql, kf = k"^ = 0, the full WW^yj vertex is thus given by 

r^^"^" = rf'"' (1 + z-fw) . (8) 

In the expressions for the two vertices, terms proportional to fcf or k2 have been dropped. 
Such terms will be contracted with the photon polarization vectors e*^{ki) or 6*'^{k2), or 
a conserved electromagnetic current and hence vanish in the amplitude. Similarly, terms 
proportional to qf are dropped since, in the massless quark limit, the W couples to a 
conserved quark current. No such assumption is made for the W-decaj leptons, and hence 
our expressions are valid when including finite charged lepton masses. For off-shell photons 



or space-like ly-bosons, more complicated expressions are obtained ^ . 

By construction, the resulting amplitude for the process qq' -^ £"1/77 should be gauge 
invariant. Indeed, gauge invariance of the full amplitude can be traced to the electromagnetic 
Ward identities 

fc^F,/ = {tDw)-l (gi) - iiDw)-^ (g2) , (9) 



between the WW'j vertex and the inverse W propagator p6[ and 

kiT^pf^uiqi, q2, ki, k2) = e (F„^^(gi, gs + ki, Ajs) - T^puiqi - ki, q2, k2)) , (10) 

relating three- and four-point functions. Since 

k,T-^^ = ({qlg-P - gf gf ) - {qlg^^ - q^q^^)) (1 + ^7w) , (H) 

and 

{'iDw)^l (g) = (g^ -ml. + iq^-fw) ( ^a/3 - ^/ J - ^w^T' ' (^2) 

the Ward identity of Eq. (^ is satisfied for the W propagator and WW'j vertex of Eqs. (||) 
and (^. Similarly the Ward identity for the WW'j'y vertex is verified for the explicit three- 
and four-point functions of Eqs. (^ and (§). 

10 



Extending this analysis to the WW''yyy vertex, the relevant Ward identity relating 
WW'j'j'j- and WW'j'j vertices is given by 

^i/^r|^/3i.p(9i5 Q2, ki, /C2, ks) = e {Tap^p{qi, q2 + ki, /c2, ks) - Tapup{qi - ki, q2, /c2, ks)) . (13) 

The right hand side vanishes for the tree-level WW'j'y vertex and thus also for the fermion- 
one-loop corrected vertex of Eq. (H). This means that the amplitude for the three photon 
process qq' -^ £^z/777 is rendered gauge invariant by implementing the corrected WWj 
and WW'yy vertex functions only, but without taking into account any WW'j'y'y one-loop 
correction. The argument can immediately be generalized to an arbitrary number of final 
state photons. For hard, non-collinear photon emission this is mostly of theoretical interest, 
however, since the cross section for 14^777 production, already, is expected to be too small 
to be observed even at a high luminosity Tevatron. 

Returning to the calculation of the qq' -^ £1/77 amplitude, a gauge invariant result 
is obtained by replacing all VT-propagators, WW'j and WW'-ff vertices in the Feynman 
graphs of Fig. |T] by the full expressions of Eqs. (^, (|^ and (^, respectively. Formally, 
these expressions include the imaginary parts of up to two loops in the vertices of Fig. |I](e). 
However, the Dyson resummation of the l^-propagators already constitutes a mixing of all 
orders of perturbation theory and thus the appearance of several vertex loops should be no 
surprise. This result is obtained naturally by attaching the two photons in all possible ways 
to either one of the fermion loops or to one of the lowest order W propagators in the zero, 
one, two etc. fermion bubble graphs contributing to the Dyson resummed process qq' -^ iu: 
the remaining sum over W vacuum polarization graphs restores the full W propagator of 
Eq. (^ on either side of a triangle graph, a WW'j vertex, the box graph, or the WW'j'y 
vertex. After resummation, one therefore obtains the Feynman graphs of Fig. |1] where any 
WW'-f{'j) vertex is given by the sum of the lowest order vertex and the imaginary part of 
the triangle (box) graphs, as defined in Figs. and ^j. 

Finally note that conservation of the final state lepton current has not been assumed 
anywhere, i.e. terms proportional to ^2 have been kept throughout. Thus our calculation 

11 



is correct for massive final state leptons and the emission of two photons coUinear with 
the charged final state lepton can be simulated with the resulting code |]12[. Alternative 
approaches in treating unstable gauge bosons in a gauge invariant way have been discussed 
in Ref . E 



III. SEARCHING FOR THE RADIATION ZERO IN W-f-f PRODUCTION AT 

THE TEVATRON 

A. Input Parameters and Detector Simulation 

We now study in detail the radiation zero in qq' — > 1^77 predicted by the SM, for pp 
collisions at ^/s = 2 TeV. To simplify the discussion, we shall concentrate on the W'^j, 
W~ —>■ e~u channel. In pp collisions, the total cross sections for W~'-f'~f and W'^'j'j pro- 
duction are equal. Angular and rapidity distributions for the W^ case can be obtained 
by a sign change of the variable. The parameters used in our numerical simulations are 
mw = 80.22 GeV, rriz = 91.187 GeV, and Oem = 1/128. We use the parton distribution 



functions set A of Martin-Roberts-Stirling |^ with the factorization scale set equal to the 
parton center of mass energy. Vs. 

To simulate the finite acceptance of detectors, we impose cuts on observable particles in 
the final state. Unless otherwise stated, we require: 

p'^ > 10 GeV, ly^l < 2.5, AR^^ > 0.3 for photons, 

Pt > 15 GeV, \ye\ < 2.5, AR^-y > 0.7 for charged leptons, (14) 

and 

^^ > 15 GeV. (15) 

Here, pt is the transverse momentum and y the rapidity of a particle, and ^rp denotes the 
missing transverse momentum of the event, defined by the imbalance to p^ and pj- iii our 
calculation. For massless particles, the rapidity and the pseudorapidity, rj, coincide. 

12 



AR = yJ{A(f)y + (Ar7)2 (16) 

denotes the separation in the pseudorapidity-azimuthal angle plane. Without finite pj^ and 
A_Re7 cuts, the cross section for ez/77 production would diverge, due to the various collinear 
and infrared singularities present. 

As mentioned in the Introduction, it is instructive to compare the results obtained for 
qq' -^ e^z/77 with those for the neutral channel, qq -^ e^e~77. In this case, we also impose 



M(e+e") > 20 GeV (17) 

cut to avoid the mass singularity from timelike virtual photon exchange graphs. The qq — »• 
e''"e~77 helicity amplitudes were calculated using the same technique which we employed in 
the 61^77 case. The transverse momentum and rapidity cuts listed above approximate the 



phase-space region which will be covered by the upgraded CDF |3J] and D0 detectors pS 

Uncertainties in the energy measurement of electrons and photons are, unless stated 
otherwise, taken into account in our numerical simulations by Gaussian smearing with 

a 0.2 

- = -^®0.01. (18) 

where the two terms are added in quadrature and E is in units of GeV. The only visible 
effect of the finite energy resolution in the figures presented below arises in regions of phase 
space where the cross section changes very rapidly, e.g. around the W or Z boson peaks. 

For the cuts listed in Eq. (0), backgrounds to ez/77 and e+e~77 production are small, 
provided the two photons are well isolated from any hadronic energy in the event. The iso- 
lation cut essentially ehminates the backgrounds from W'-y + 1 jet and ly + 2 jet production 
where one or both jets fragment into a photon [^. For pj^ > 10 GeV, the probability that 
a jet fakes a photon, Pj/-y, is 10^^ or less 0]. Backgrounds from W'y+ jets and W+ jets 
production, where one or two jets fake a photon, are then small. The photon-photon sepa- 
ration cut of AR^^ > 0.3, combined with a substantial py, requires a sizable invariant mass 
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of the two-photon system and thereby ehminates backgrounds from 7r° -^ 77 decays which 
might originate from W/Z + 1 jet production with a leading vr". 

The geometrical acceptance of the upgraded CDF and D0 detectors for muons will be 
similar to that for electrons. Requiring the charged lepton to be well separated from the 
photons, the cross sections for 61/77 and fiw-f-f production are then nearly identical. The 
results derived in the following for the electron channel therefore also apply to the fiwyy 
final state. 

B. The e 1^77 and W^{—>- e^u)^^ Cross Sections 

In Fig. ^, we present the total cross sections, within the cuts of Sec. IIIA, for pp — > e~z/77 
and on-shell W~{-^ e~z/)77 production (solid) as a function of the pp center of mass energy. 
For comparison, we also show the cross sections for e"'"e~77 and Z{-^ 6+6^)77 production 
(dashed). Here the on-shell W~{~^ e~u)'~f'~f and Z{—>- 6^6^)77 cross sections have been 
calculated in the narrow W/Z width approximations. The large differences between the on- 
shell and the full cross sections arise from diagrams where one or both photons are emitted 
by a final state charged lepton. For e"^e^77 events there are also sizable contributions from 
7* -^ e^e" . Contributions from these diagrams increase the cross section by about a factor 
3 (6) in the 14^77 (^'77) case for the cuts chosen. No energy smearing effects are taken into 
account in Fig. |[ 



The rates for pp -^ W~{-^ e~z/)77 and pp -^ Z{-^ e"'"e~)77 [Tj] are almost identical 
over the entire center of mass range studied for the cuts chosen. This should be contrasted 
with the cross section ratio of W~{—* e~z/) + 2 jet to Z{-^ e^e") + 2 jet production which 



is about 4.6 ||3^. The relative suppression of the ^^77 cross section can be traced to the 
radiation zero which is present in Wj'-f, but not in Z77 production. Similarly, the W'~f cross 
section is suppressed relative to the Z'j production rate because of the radiation zero in 
qq'^Wjm- 



Figure || shows that, although we require the charged lepton to be well separated from 
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FIG. 4. The total cross sections for pp -^ e~u^^ and pp -^ W^{-^ e^ 1^)77 (solid lines) as 
a function of the pp center of mass energy, ^/s. For comparison, we also show the pp -^ e^ e~^^ 
and pp -^ Z(-^ e^e~)77 cross sections (dashed lines). The acceptance cuts are summarized in 
Sec. IIIA. No energy smearing is imposed. 
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the photons, radiation off the final state charged lepton completely dominates. In order 
to search for a possible radiation zero in W'-ff production, it is necessary to suppress final 
state bremsstrahlung more efficiently. To isolate the W{—>- eu)'-/'-/ component in pp -^ ewyy, 
it is useful to study the transverse mass distribution of the eu system which is shown in 
Fig. §(a). W{^ ei/)77 events produce a Mxieu) distribution which is sharply peaked at 
Mr(ei/) = rriw- However, finite energy resolution effects significantly dilute this peak (see 
solid curve). On the other hand, if one or both photons are emitted by the charged lepton, 
the eu transverse mass tends to be considerably smaller than the W mass. Requiring 

MT{e~iy) > 70 GeV, (19) 

eliminates most of the contributions from final state radiation. With this additional cut, 
the total and differential cross sections for 61/77 and W{-^- 61^)77 production are almost 
identical. 

Similarly, a cut on the di-lepton invariant mass can be used to suppress photon radiation 
from the final state leptons in pp -^ 6^6^77. The e"''e~ invariant mass distribution, for the 
cuts summarized in Sec. IIIA, is shown in Fig. ^(b). The two broad peaks below the Z 
resonance region correspond to contributions from Z -^ e^e~7 and Z -^ e"'"e~77. Details 
of the structure depend on the choices of pj^ and ARe^ cuts. For 

M(e+e") > 85 GeV, (20) 

contributions from final state bremsstrahlung are reduced by about a factor of 4 for pj, > 
10 GeV. Nevertheless, contributions from final state bremsstrahlung and 7* — > e^e~ are still 
sizeable in this case. 

Within the cuts of Eqs. (|1^) and (p!5|), and with the transverse mass cut of Mt{(!-v) > 
70 GeV, the total ^^z/77 (£ = e, /i) cross section for pp collisions at ^/s = 2 TeV is about 
2 fb. For an integrated luminosity of 30 fb~^, one thus expects about 60 £^1/77 events. 
The cross section depends quite sensitively on the minimum photon transverse momentum, 
j9™™, however. This dependence, with and without the transverse mass cut of Eq. ([T9|), is 
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FIG. 5. (a) The ev transverse mass distribution for pp — > e~i>77, and (b) the e~^e~ invariant 
mass distribution for pp -^ e^e~77, at ^/s = 2 TeV. The cuts imposed are summarized in Sec. IIIA. 
The solid and dotted curves give the results with and without taking into account the finite energy 
resolution of detectors (see Eq. (p^)). 
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shown in Fig. ^ for the e~i>77 cross section. For completeness, curves for pp -^ e^e~77 
are included as well. No energy smearing effects are taken into account here. Reducing the 
photon transverse momentum threshold from 10 GeV to 4 GeV, the e~i^77 rate, regardless 
of the transverse mass cut of Eq. (|T9|), increases by about a factor of 6. Due to the limited 
number of ez/77 events even at the highest Tevatron luminosities, the pt threshold for at 
least one of the photons should be lowered as far as possible in a search for the radiation 
zero in 1^77 production. Nevertheless, in our further analysis, we shall retain the more 
stringent photon transverse momentum requirement of pj^ > 10 GeV for both photons (see 
Eq. (p!^) ). As mentioned in Sec. IIIA, backgrounds from 14^7+ jets and W+ jets production, 
where one or two jets fake a photon, are then small. Furthermore, we shall impose the mass 
cuts of Eq. (|19D and pO|) unless stated otherwise. 

C. Searching for the Radiation Zero 

The general theorem of Ref. states that, in the SM, the amplitude for the process 
du -^ iy~77 vanishes for 

cos^^^ = cos^ot^- = _ = --, (21) 

when the two photons are collinear. Here 6'^ is the angle between the incoming d-quark 
and the W boson, and the asterisk on a quantity denotes that it is to be taken in the parton 
center of mass frame. For W^'^'^ production, the role of the u- and rf-quarks in Eq. (|2lD 
are interchanged, i. e. cos 0^^^+ = — cos 9Q^r- . The existence of the radiation zero can be 
readily verified numerically. Figure ^(a) shows the cos 6^ distribution for the parton level 
process du -^ W'^yy, at a parton center of mass energy of V^ = 300 GeV and for M^^ = 0, 
which forces the two photons to be collinear. In addition, the photon energies are chosen to 
be equal. For unequal photon energies, qualitatively very similar results are obtained. The 
vanishing of the differential cross section at cos6'Qy(^_ = —1/3 is apparent. For comparison, 
we have also included the cos 6*^ distribution for uu — > Z77 in Fig. |^(a) (dashed line), and 
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FIG. 6. The total cross sections for pp -^ e~v^^ and pp -^ 6^6^77 at y/s = 2 TeV as a function 
of the minimal photon transverse momentum p™™. The solid lines are for the cuts summarized in 
Sec. IIIA. For the dashed lines an additional Mxieu) > 70 GeV or M(e'''e~) > 85 GeV cut has 
been imposed. No energy smearing is included here. 
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FIG. 7. The angular distribution of the vector boson for the partonic processes (a) 

du -^ VF~77 with Myy = 0, and (b) du -^ W^'j, at \/s = 300 GeV (solid lines). Corresponding 
curves for uu -^ Zjj and uu -^ Z7 (dashed lines) are also shown for comparison. The photon 
energies in the 1^77 and 2^77 case are chosen to be equal. The differential cross sections are in 
arbitrary units. No cuts and no energy smearing are imposed, and the W and Z bosons are treated 
as stable particles. 
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show the cos 9^ and cos 9*^ distributions for du — > W~'y and uu -^ Z'j at yl = 300 GeV in 
Fig. |^(b). The cos^^ region where the cross section is substantially reduced due to the zero 
is seen to be considerably larger in the W''yy case. No cuts and no energy smearing have 
been imposed in Fig. |^, and the W and Z bosons are treated as stable particles. The overall 
normalization of the cross sections in each part of the figure is arbitrary. Similar results are 
obtained for different parton center of mass energies. 

The impressive W'j'j radiation zero in Fig. |^ is washed out by the small contamination 
of W{^ e~z/7)7 and W -^ e~u'yy events which pass the Mxie^i^) cut of Eq. (|19|) when W 
decays and finite W width effects are taken into account. Binning effects reduce the radiation 
zero to a mere dip as well. This is shown in Fig. ^ where we display the normalized double 
differential cross section (l/o") ((i^o"/(i cos 6'^ c? cos 6** ) for the partonic process du —>■ e'u'yy 
with vS = 300 GeV. Here 9* is the angle between the two photons in the parton center 
of mass frame. In this figure, the full set of contributing Feynman diagrams, including the 
corrections to the W propagator and the WW'j and WWj'j vertices described in Sec. II, 
have been taken into account, together with the cuts summarized in Eqs. (p!^), (|I5[) and (|19|) 
which we subsequently impose in all figures. 

Figure || demonstrates that the dip in the e~z/77 differential cross section at cos 9^ = 
— 1/3, which signals the presence of the 1^77 radiation zero, is quite pronounced for the 
cuts we have chosen. It also shows that it only gradually vanishes for increasing values 
of 9* . Requiring two photons with AR^^ > 0.3 therefore has no significant effect on the 
observability of the W'j'j radiation zero. 

The significance of the dip, which signals the amplitude zero, is potentially further re- 
duced by the convolution with parton distribution functions and by the twofold ambiguity 
in reconstructing the parton center of mass frame. This twofold ambiguity originates from 
the non-observation of the neutrino arising from W decay. Setting the eu invariant mass 
equal to mw leaves two solutions for the reconstructed center of mass, which can be ordered 
according to whether the rapidity of the neutrino is larger ( "plus" solution) or smaller ( "mi- 
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FIG. 8. The normalized double differential distribution (l/cr) [cP a / d cos 9* d cos 6^^) for the 



partonic process du ^ e z/77 at v^ = 300 GeV. The cuts summarized in Eqs. (p^, (|T5|) and (19) 
are imposed. 
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nus" solution) than the rapidity of the electron . Since the photons couple more strongly 
to the incoming up-type anti-quark, the W~ boson tends to be emitted in the proton direc- 
tion. Within the SM, and as in W'j production, the dominant helicity of the W~ boson in 
W'-f'j production is Xw = — 1, implying that the electron is more likely to be emitted in the 
direction of the parent W. The rapidity of the electron thus is typically larger than that 
of the neutrino, and the "minus" solution better preserves the dip caused by the radiation 
zero. In W'^'yy production, the W boson is dominantly emitted into the p direction, and, 
consequently, the "plus" solution shows more similarity with the true reconstructed parton 
center of mass. 

The normalized double differential distribution (1/cr) {d'^a/d cos 9*^ dcos9^) for the pro- 
cess pp -^ e~iJ'~f'~f at y/s = 2 TeV is shown in Fig. ^, using the "minus" solution for the 
reconstructed parton center of mass. The distribution is seen to be quite similar to the 
corresponding partonic differential cross section shown in Fig. §. The convolution with the 
parton distribution functions therefore has only a minor effect on the observability of the ra- 
diation zero. Likewise, the reconstruction of the parton center of mass frame does not affect 
the significance of the dip much, provided that the appropriate solution for the longitudinal 
momentum of the neutrino is used, and the missing transverse momentum is well measured 
(see below). 

For the limited number of ewyj events expected in future Tevatron runs, it will be 
impossible to map out the double differential distribution shown in Fig. ^. However, since 
the dip signaling the radiation zero disappears only gradually with increasing values of 
9*^, most of the information present in (i^o"/(i cos 6'*^ (i cos 6*^ is contained in the cos^^ 
distributions for events with cos^*^ > versus cos^*^ < 0. These two cos^^ distributions 
are shown in Fig. 10, for both the "plus" and the "minus" solution of the reconstructed 



parton center of mass frame. For comparison. Fig. |T0] also shows the cos 9*^ distribution 
for pp — > e"'"e~77 with M{e~^e~) > 85 GeV and the cuts of Eq. ([I^). For cos 6**^ > and 
using the "minus" solution, the cos 6*^ distribution displays a pronounced dip located at 
cos^^ ~ —1/3. For the "plus" solution the minimum is shifted to cos9^ ~ 0. In contrast, 
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FIG. 9. The normalized double differential distribution (l/cr) {d? a / d cos 9* d cos 9^r) for the 
process pp — > e'uj'j at ^/s = 2 TeV, using the "minus" solution for the reconstructed parton center 
of mass. The cuts summarized in Eqs. (p!4|), ( |l5|) and (19) are imposed. 
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comparison. The cuts summarized in Eqs. (14), ( p^ ) and (|T9|) are imposed. In the .2^77 case, the 
cuts listed in Eqs. (Ill) and (pH) are applied. 



requiring cos^* < 0, the dip is drastically reduced, and the differential cross section at 
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cos 6'^ ~ —1/3 is about one order of magnitude larger than for cos6'*^ > (see Fig. |TO|(b)). 
The large difference in the cos 9^ distribution for cos 9*^ > and cos 9*^ < becomes more 
apparent by comparing the cos 9^ distribution in Z77 production in the two regions. Unlike 
the situation encountered in W'jj production, the cos^^ distributions for cos 6**.^ > and 
cos 6'*.^ < are very similar. 

Determining the cos 9^ distribution requires measurement of the transverse momentum 
of the neutrino produced in the W decay. In 61^77 production, the neutrino transverse mo- 
mentum is identified with the missing transverse energy, ^^, in the event. In future Tevatron 
runs, one expects up to ten interactions per bunch crossing ^^. Multiple interactions per 



25 



crossing significantly worsen the Ifj, resolution, and thus tend to wash out the dip signaling 
the radiation zero. We have not included missing transverse energy resolution effects in our 
simulations, since the number of interactions per crossing, and hence the Ifrp resolution, sen- 
sitively depend on the future Tevatron accelerator parameters which are difficult to foresee 
at present. 

Due to the negative impact of multiple interactions on the missing transverse energy 
resolution, it is advantageous to search for a kinematic variable which exhibits a clear signal 
of the radiation zero but does not depend on the neutrino momentum. The y^^ — ye distri- 
bution is a possible candidate for such a variable. Here ye is the electron rapidity and y^^ 
denotes the rapidity of the two-photon system in the laboratory frame. 

In Ref. I^H it was found that photon lepton rapidity correlations are a useful tool to 
search for the radiation zero in W'-f production. The distribution of the rapidity difference, 
A|/(7, e) = y.y — ye, which does not require knowledge of the missing transverse energy or 
the longitudinal momentum of the neutrino, clearly displays the SM radiation zero in form 
of a dip. In the parton center of mass frame, the photon and W boson in qq' -^ W~'y are 
back to back. Due to the radiation zero, the photon and W rapidity distributions in the 
parton center of mass frame, da/dy* and da/dy^, display pronounced dips located at 

y; = ^ln2^0.35, (22) 

y*w - -0-05. (23) 

If W mass effects could be ignored, one would expect that y* = —y'w- Since differences of 
rapidities are invariant under longitudinal boosts, the difference of the photon and the W 
rapidity in the laboratory frame then exhibits a dip at 

A|/(7, W) = y,-yw = y;-y*w- 0.4. (24) 

As discussed earlier, the dominant W helicity in iy^7 production is \w = ±1, implying 
that the charged lepton tends to be emitted in the direction of the parent W , and thus 
reflects most of its kinematic properties. The dip signaling the presence of the radiation zero 
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therefore manifests itself in the A|/(7, £) distribution. Since the average rapidity of the lepton 
and the W are shghtly different, the location of the minimum is shifted to A?/(7, £) ^ 0.1. 

The radiation zero in W'-ff production occurs at exactly the same rapidity as the zero in 
Wj production, when the photons are coUinear. One therefore expects that the Ay(77, W) 
distribution displays a clear dip for photons with a small opening angle, 6^-y, in the labo- 
ratory frame, i.e. at cos^^^ ^ 1. In Fig. |ll| we show the double differential distribution 
d?a /dcosO^^ d{yyy — yw), using the "minus" solution for the longitudinal neutrino momen- 
tum. For A|/(77, 14^) ^ 0.4, a clear dip is visible for cos 6*^^ values close to one. The dip 
gradually vanishes for larger opening angles between the two photons, leading to a "canyon" 
in the double differential distribution. Due to the finite invariant mass of the 77 system 
for non-zero values of 6^^, the location of the minimum in A|/(77, W) varies slightly with 
cos 9^^. 

Since the dip vanishes gradually with decreasing cos 6'^^, it is useful to consider the 
Ay(77, W) distribution for cos^^^ > and cos^^^ < 0. Figure |12](a) displays a pronounced 
dip in da/dAyi^'yy, W) for cos^^^ > 0, located at Ay (77, W) ~ 0.7 (solid line). In contrast, 
for cos 6*^^ < 0, the Ay{'yj,W) distribution does not exhibit a dip (dashed line). The 
A|/(77, W) distribution for cos 6*^^ > thus plays a role similar to the Ay(7, W) distribution 
in W'j production. In the dip region, the differential cross section for cos^^^ < is about 
one order of magnitude larger than for cos 6'^^ > 0. In addition, the Aj/(77, W) distribution 
extends to significantly higher y^^ — yw values if one requires cos6'^^ > 0. This reflects the 
narrower rapidity distribution of the two-photon system for cos^-y^ < 0, due to the larger 
invariant mass of the system when the two photons are well separated. 

Exactly as in the W'y case, the dominant helicity of the W boson in 14^^77 production 
is Xw = ±1- One therefore expects that the distribution of the rapidity difference of the 
77 system and the electron is very similar to the y^^ — yw distribution and shows a clear 
signal of the radiation zero for positive values of cos6'^-y. The y^^ — ye distribution, shown in 
Fig. |T^(b), indeed clearly displays these features. Due to the finite difference between the 
electron and the W rapidities, the location of the minimum is again slightly shifted. The 
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FIG. 11. The double differential distribution d?a/dcos9^^d{y^^ — yw) for the process 
pp — > e~P77 at ^/s = 2 TeV, using the "minus" solution for the longitudinal neutrino momen- 
tum. The cuts summarized in Eqs. (0), (^) and (|l9|) are imposed. 
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Ai?^^ > 0.3 cut has only little effect on the significance of the dip. 

The characteristic differences between the A?/(77, e) = y^^—i/e distribution for cos 6'^^ > 
and cos^^^ < are also reflected in the cross section ratio 

n = {^^(^^'-)>-^ f , (25) 

which may be useful for small event samples. Many experimental uncertainties cancel in TZ. 
For cos6'^^ > one finds TZ ~ 0.25, whereas for cos6'^^ < 7?. ~ 1.06. 

Our calculations have all been carried out in the Born approximation. The complete 
NLO QCD corrections to W'-ff production have not been calculated yet; only the hard 
jet corrections to W{^ ez/)77 production are known [^]. It is reasonable, however, to 
take the known NLO QCD correction to W'y production as a guide [|l^,|l^]. At 0{as) the 
virtual corrections only enter via their interference with the Born amplitude, and thus the 
radiation zero is preserved in the product. Among the real emission corrections, quark- 
antiquark annihilation processes dominate at Tevatron energies. According to the theorem 
of Ref. 1^, extra gluon emission, i.e. the process qq' -^ W'^n'yg, exhibits a radiation zero 
at cos 6*^ = ±1/3 if the gluon is collinear to all emitted photons, and also in the soft gluon 
limit. Eg ^ (again, provided the photons are collinear). This leaves quark-gluon initiated 
processes to potentially spoil the radiation zero. They are still suppressed at the Tevatron, 
however, especially when a large photon-jet separation is required. As a result, we expect 
the dip signaling the radiation amplitude zero to remain observable, at Tevatron energies, 
once NLO corrections are included. 

At the LHC {pp collisions at ^/s = 14 TeV [0), the bulk of the QCD corrections to 
W'j'y production originates from quark gluon fusion and the kinematical region where the 
final state quark radiates a soft W boson which is almost collinear to the quark. Events 
which originate from this phase space region usually contain a high pt jet. Since there is 
no radiation zero present in the dominating qg — > W'-ffq' and gq' —>■ W'-ffq processes, it is 
likely that QCD corrections considerably obscure the signal of the W'-ff radiation zero at the 
LHC, as in the Wj case [^. This conjecture is supported by the large relative cross section 
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of VF77 + 1 jet production as compared to W'-yy production reported in Ref. p9[. Although 
a jet veto should help reducing the size of the QCD corrections, NLO QCD corrections to 
W^77 + jet production may still significantly reduce the observability of the radiation zero 
for jet definition criteria which are realistic at LHC energies. We therefore do not consider 
62^77 production at the LHC in more detail here. 

IV. SUMMARY AND CONCLUSIONS 

We have presented a calculation of the process pp -^ ez/77 including final state 
bremsstrahlung diagrams and finite W width effects, and explored the prospects to ob- 
serve the radiation zero predicted by the SM for pp -^ Wyj in future Tevatron collider 
experiments. In order to obtain a gauge invariant scattering amplitude, the imaginary parts 
of the WlVj triangle graphs and WW'-ff box diagrams have to be included, in addition 
to resumming the imaginary contributions to the W vacuum polarization. The imaginary 
parts of the triangle and box diagrams were found to change the lowest order WWy and 
WW'-ff vertex functions by a factor [l + iTw/rnw) for the momentum configuration relevant 
for the process qq' -^ ew-f-f. A gauge invariant result for the qq' -^ ez/77 amplitude is then 
obtained by replacing all W propagators, WW'-y and WW'-yy vertices by the full expressions 
of Eqs. (H), (0) and (§), respectively. The same prescription also ensures that the Ward 
identities relating the WWnj and WW{n — 1)7, n > 3, vertex functions are fulfilled, and 
thus yield a gauge invariant amplitude for qq' —>■ eu + nj with n > 3, without taking into 
account one-loop corrections to these higher vertex functions. 

The SM predicts the existence of a radiation zero in qq' -^ W^'-yy at cos6'^ = ±1/3 if 
the two photons are collinear. Here 6^ is the angle between the W and the incoming quark 
in the parton center of mass frame. Since it is very difficult to experimentally separate 
two collinear photons, one has to search for a signal of the radiation zero which survives 
an explicit photon-photon separation requirement. Contributions from Feynman diagrams 
where one or both photons are emitted by the final state charged lepton eliminate the 
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radiation zero and therefore need to be suppressed by suitable cuts. We found that a large 
lepton-photon separation of A/?e^ > 0.7, together with a cut on the ev transverse mass of 
Mr(ei/) > 70 GeV suppresses these contributions sufficiently. 

The 14^77 radiation zero is signaled by a pronounced dip in the cos^J^ distribution 
if one requires cos^* > 0. In contrast, no dip is present for cos0* < 0. In order to 
measure the cos ^^ distribution, the parton center of mass frame has to be reconstructed. 
Since the neutrino originating from the W decay is not observed in the detector, this is 
only possible modulo a twofold ambiguity. The two solutions can be ordered according to 
whether the reconstructed rapidity of the neutrino is larger ("plus" solution) or smaller 
("minus" solution) than the rapidity of the charged lepton. For iy~77 {W^'^'^) production, 
the "minus" ( "plus" ) solution is found to best represent the expected kinematical features. 

When searching for the radiation zero in W^77 production it is advantageous to consider 
alternate variables which, unlike the cos^^ distribution, do not depend on the neutrino 
momentum. The rapidity difference between the two-photon system and the electron, i.e. 
the y-y^ — He distribution, fulfills this requirement. It was found to exhibit a pronounced dip 
which signals the presence of the radiation zero if a cos 6^^ > cut is imposed {6^^ being the 
opening angle between the two photons in the laboratory system) . As expected, the y^^ — ye 
distribution shows no dip for cos 6*^^ < 0. A photon-photon separation cut of AR^^ > 0.3 
has little effect on the observability of the radiation zero. Although we have restricted our 
discussion to ewyy production, our results also apply to pp — > fiwyy. 

The conditions for which one expects a radiation zero in the SM qq' -^ W''yy and qq' — >■ 
W'j amplitudes and the location of the zeros are closely related: the four-momentum of the 
photon in W'y production simply has to be replaced by the four-momentum of the 77 system 
in the W'j'-f case with the additional requirement that the two photons are coUinear. We 
have demonstrated that a similar replacement in the W'j photon-lepton rapidity difference 
distribution, with the less stringent requirement on the opening angle between the photons 
of cos 6*^^ > 0, is in fact sufficient to produce an observable signal of the 1^77 radiation zero 
(see Fig. |l|). 
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NLO QCD corrections to pp -^ W'-ff are expected to be modest at Tevatron energies. 
Given a sufficiently large integrated luminosity, experiments at the Tevatron studying cor- 
relations between the rapidity of the photon pair and the charged lepton therefore offer an 
excellent opportunity to search for the SM radiation zero in hadronic W'-ff production. 
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